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Abstract 



The p, K* and <j) decay constants and the vector current renormalization 
factor are studied by using an 0(a 2 ) classically-improved, tadpole-improved 
action. Tree-level calculations are used to show how the classical improvement 
of the action, involving next-nearest-neighbour timesteps, is transferred to the 
matrix elements. Simulations are performed on coarse lattices and compared 
to Wilson results from both coarse and fine lattices. The improved action 
data are found to resemble Wilson data obtained at 1/3 of the lattice spacing, 
which is the same degree of improvement that is seen by comparing the mass 
spectra. 
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I. INTRODUCTION 



Numerical simulations of QCD can be performed by discretizing Euclidean space-time 
onto a lattice and working within a finite volume. Physical measurements are then obtained 
by extrapolating to zero lattice spacing and infinite lattice volume. There are many lattice 
actions which correspond to continuum QCD in the appropriate limit, and a lot of recent 
work has focussed on finding actions which reduce the errors due to discrete lattice spacing. 
The goal is to find actions which provide accurate results on coarsely-spaced lattices and so 
allow lattice sizes (and therefore computing resources) to be reduced. A progress report on 
improving actions may be found in Ref. [|I|]. 

To demonstrate the value of an improved action, one would like to show that accurate 
results can be obtained from coarse lattices for a wide variety of quantities. As a first step 
beyond the mass spectrum, where significant improvements on coarse lattices have already 
been demonstrated, we consider the decay constants of light vector mesons and the associated 
vector current renormalization factor. (A preliminary discussion was given in Ref. [0].) An 
advantage of studying vector currents is the existence of a conserved vector current even at 
nonzero lattice spacing, which does not get renormalized. Neglecting lattice spacing errors, 
the matrix element of the conserved current between any initial and final states is identical 
to the matrix element of the local vector current between those same states, except for the 
multiplicative "renormalization factor", Zy. In this way, Zy can be easily computed. In 
practice, it has been known for a long time || that the determination of Zy is quite sensitive 
to lattice spacing errors, which is another reason that the vector current is valuable for the 
present work, since it is precisely these errors that should be reduced by an improved action. 

The action [|],[5| that will be used here is improved at the classical level up to and including 
0(a 2 ), and is also improved at the quantum level via the "tadpole factor" of Lepage and 
Mackenzie ||. The action has the added feature of being quite simple, since the gauge 
sector contains only lxl and 1x2 plaquettes, and the fermion sector does not contain the 
so-called "clover" term and has no need of the associated field transformations. 

The mass spectrum of this improved action has been studied numerically in Ref. J?|, 
where it was concluded that data from the original Wilson action || at a lattice spacing 
of a ~ O.lfm can be reproduced with the improved action on lattices with a ~ 0.3fm, 
representing a very substantial saving of computing resources. Other studies of the mass 
spectrum with the same improved action have produced similar results |]9|JT0|1. 

In this paper it is shown that a comparable degree of improvement is also obtained in the 
computation of the vector meson decay constants, f p , fx* and f^, and the renormalization 
factor, Zy. Sec. [Ij begins with a review of the general methods for calculating vector 
current matrix elements, and then discusses the particular improved action chosen for this 
work. Sec. |T| is used to show explicitly that the quantities to be computed are classically 
improved. The results of numerical simulations using the improved action are compared to 
Wilson action results in Sec. |IV|, and a summary is given in Sec. [V]. 
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II. METHOD OF CALCULATION 



The p meson decay constant, f p , is denned from a continuum matrix element of the local 
vector current as follows. 
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(0|F» cont = (1) 
Jp 



V^{x)=^{x) 1 ^{x) (2) 

A similar definition applies to fx* and f,p. Recall that this current is conserved in the 
continuum. 

At nonzero lattice spacing the local current is not conserved, and its matrix element 
between any initial and final states gets renormalized by a factor Zy, 

(f\V^) cont = Z v (g 2 )(f\V^) + 0(o n ) • (3) 

g is the lattice coupling, and 0(a n ) represents the contributions from operators of higher 
mass dimension which are multiplied by powers of the lattice spacing, a. For the Wilson 
action these corrections begin at 0(a), but for our improved action no 0(a,a 2 ) classical 
terms can appear, and the tadpole factor is expected to approximately minimize terms 
proportional to the coupling, for example 0(g 2 a) terms. 

If all quark flavours have the same mass, then both the improved action and the Wilson 
action have an exact vector (flavour) symmetry, which implies the existence of a conserved 
"Noether" current, V?f . (Actually, there is not a unique conserved current, and care must 
be used to choose one that is classically improved to the appropriate order. This will be 
addressed below.) Because of the nonrenormalization of conserved currents, 

(fK L ^)cont = (f\V^) + 0(a n ) , (4) 

which allows Zy to be computed from a ratio of matrix elements as follows. 

Zv{g2) = mm + ° {an) (5) 

Although Zy itself is independent of the initial and final states, the ratio of matrix elements 
in Eq. (|5|) is not, so the size of the lattice spacing errors can be estimated by varying \i) and 
\f). In this work, we will consider four separate choices for the ratio in Eq. (^): a vector 
current insertion between pseudoscalar meson, vector meson or spin-1/2 baryon states, 

R _ (PW°\P) R _ (VK°W) R _ (B\Vf\B) 

p ~(pK\p) ' v ~(y\v, L \v) ' Kb ~ (B\Vt \B) w 

and the decay process, 

B _ (QKW) m 

Rd = WW)- () 
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In our lattice simulations, states will be created via the following local interpolating 
fields. 



pseudoscalar mesons : x ( x ) = i^i^l^i^) 
vector mesons : xZ{ x ) = 4'( x )'l^'4 ) ( x ) 



spin—- baryons : X B ( X ) 



(8) 
(9) 

(10) 



(C is the charge conjugation matrix.) The decay constant itself is obtained from simulations 
of the two-point vector correlator, 



^m 3 f>~ m vt 



for t > 0, 



(11) 



which is derived by multiplying Eq. ([I]) with its hermitian conjugate, converting from con- 
tinuum to lattice quantities, and performing the appropriate summations. The explicit de- 
pendence on Zy can be changed by replacing one (or both) of the local currents in Eq. ( |TT|) 
by a conserved current. On a lattice with periodic boundary conditions, another term must 
be added to the right-hand side of Eq. (|TT|) which decays exponentially from the (periodic) 
image of the source. 

Following the work of Luscher and Weisz 0], the gauge field term of the improved action 
involves a sum over 1x2 rectangular plaquettes (U rt ) as well as 1 x 1 elementary plaquettes 
(U P i), 



S G (U) = ^ReTr 



E(i 

. p 1 



i 



pi) 



20t/ 2 rt 



E(l " Urt) 



(12) 



where a tadpole factor, defined by 



U = (^ReTrf^) 1 / 4 



(13) 



has been introduced to absorb the lattice tadpole effects and thereby improve the matching 
to perturbation theory ||. The fermion part, which was first introduced more than a decade 
ago ||, also contains next-nearest-neighbour interactions. 



X 



+ o« E [V'OXl - ln)U^x)^{x + n) + ${x + //)(! + ^)Ul(x)ij(x) 



x,ii 



6U x ^ 



Y, [$i x ){2 ~ l^U^U^x + n)ij)(x + 2/i) 



+^{x + 2/i)(2 + llt )UUx + fi)UUx)^(x) 



(14) 



The presence of next-nearest-neighbour timesteps is a generic consequence of classical im- 
provement beyond 0(a), and generates well-known artifacts, e.g. unphysical branches in 
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the free dispersion relations, discussed for gauge fields in Ref. |13[]. The effect is also seen 
for fermions, as discussed in following sections. 

The Noether vector current corresponding to Eq. (Ill) is 

Vf(x) = -$(x + /i)(l + ltx )Ul(x)4>(x) - yp{x){\ - lp )U^x)^{x + n) 



+ 



12U 
1 



+ 



12U 
1 

12U 
1 

mTn 



$(x)(2 - lp )U p {x)U p {x + n)4>(x + 2/x) 
ip(x - fi)(2 - j^U^x - fijUnix^x + n) 
${x + 2/i)(2 + l^U^x + ii)Ul{x)^)[x) 
+ fj)(2 + -yM(x)Ul(x - /i)ip(x - n) 



(15) 



Although this current is conserved, an explicit calculation shows that it is not classically 
improved to the level of the action because of extra total derivative terms. An improved 



conserved current could be written as follows UTT 
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V^Hx) = —V^(x) + —V°(x + —V°(x + u) 
13 26 26 



(16) 



but the offending total derivatives in Vf^ do not contribute to the quantities we wish to 



compute, as will be shown explicitly in Sec. |T|, so is not required for this work. 



III. CONFIRMING TREE-LEVEL IMPROVEMENT 



The matrix elements that are used to obtain f p , fx* 



and Zy can be evaluated 



analytically at the classical level. We will follow closely the approach and notation of El- 
Khadra, Kronfeld and Mackenzie ||12|| , paying particular attention to the effects of next- 



nearest-neighbour interactions, which were not present in the actions of Ref. 12 



A. The quark equation of motion 

The equation of motion that comes from Eq. ([14]) can be written in the form 



E 



527pKp(z,y) + L(x,y) 



(17) 



where the Fourier transforms of A" M (x, y) and L(x, y) are 



iKpip) = -sinp^a (4 - cosp^a) 

4 1 

L(p) = Ma + 4 - - cosp^a + - ^ cos2p M a 



(18) 
(19) 



5 



respectively. The on-shell condition is 



(20) 



which is a quartic equation in coshEa = cosj^a. Fig. |T] shows the solutions for 3-momenta 
in the (1,1,0) direction and M — 0, compared to the Wilson action result and the continuum 
E = |p] line. Notice that the improved action has artifacts which are restricted to the region 
of large Ea, and that the low-energy solution is closer to the continuum for a greater range 
of |p|a than is the Wilson curve. 

Also related to our discussion of decay constants is the limit of vanishing 3-momentum, 
p = 0, and Fig. |2] shows the solutions of Eq. ( |2"0"D versus Ma. Again, improvement is found 
with respect to the Wilson action. In this well there are next-nearest-neiehbour 

timestep artifacts which, however, remain at high energies for the values of Ma that we will 
use in numerical simulations. Note that Eq. ( p0|) is only quadratic in the mass, and becomes 
the following unique expression for Ma when p = 0. 



Ma + 1 



1 + 4cosh.Ea — 2cosh 2 Ea + (4 — coshi?a)sinh|£ , a| 



We now seek spinor solutions to the equations of motion, 



£ 7 ^(p)+£(p) 

-^7A(P)+^) 



«(£,p) = o 
«(£,p) = o 



(21) 

(22) 
(23) 



where £ = 1, 2 is the spin index. The equation of motion for u-spinors has been modified so 
that from now on E is always a positive quantity. The solutions are 



where 



"(1,0) 



w(f,p) 



A\ 
o 



voy 



^2L(p) [L(p) + sinhEa (4 - coshEa) /3] 

ilpKpjp) + L(p) 

2L(p) [L{p) + sinhEa (4 - coshEa) /3] 



(24) 
(25) 



u(2,0) 



1 




v(hO) 





1 

voy 



v(2,0) 






vi y 



(26) 



The general solution of Eq. (|17D is 



(27) 
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FIG. 1. The energy-momentum relation for a massless, free fermion from the improved action 
with 3-momentum p in the (1,1,0) direction. Only the real parts of complex solutions are shown. 
The dashed(dotted) curve is the continuum(Wilson) result. 
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FIG. 2. The energy-mass relation for a stationary, free fermion from the improved action. Only 
the real parts of complex solutions are shown. The dashed(dotted) curve is the continuum (Wilson) 
result. 
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where &(£, p) annihilates a quark and p) creates an antiquark. The normalization factor 
N(p) can be determined, for example, from the normalization of the quark propagator. By 
choosing the creation and annihilation operators to be normalized in the following way, 

{&(£', pV^Cp)} = HCPO^CP)} = (2tt) 3 5 3 ( P ' - p)#' , (28) 
the normalization factor is found to be 



(4cosLEa — cosh2i?a) sinLEa (4 — cosh^a) + 6 (2smhEa — sinh2i?a) L(p) 

(29) 



B. fv and Zy 



The vector meson decay amplitude and the insertion of a vector current between hadron 
states can now be calculated at the classical level using the matrix elements 

(OlWCpMCp')) and (?(£', p') | W£, p)) 

respectively. In each case, kinematics are chosen that corresponds to the numerical sim- 
ulations of Sec. |TV| , which means that the decay amplitude is studied in the rest frame 
(p' = — p), whereas the vector current insertion has no 4- momentum transfer (p' = p). A 
consequence of these choices is that only fi = 1,2,3 are nonzero for the decay amplitude, 
and only jj, = 4 is nonzero for the vector current insertion. 
Using as the vector current, the decay amplitude is 



(0|^%(£,p)?(£',-p)> 
= iv 2 ( P M£',-p) 

= iv 2 ( P r4 

= «(£,0) 7j u(£,0 



-cospjd — -cos2pja 



0(a 3 



4z 4 2i 1 

-—smpjd + -jjcospjd + —sm2pja — -jjcoszpja 

1 



u(e,o) 7i «(e,o) 



w (e,p) (30) 

(31) 
(32) 



where Eq. (|21| ) has been used. A similar derivation shows that Eq. (|32|) also holds for the 
local vector current, V } L , so the decay amplitude is classically improved up to 0(a 2 ) for both 
vector currents. 

The expression for the current insertion is 



<^',pM c |g(£,p)> 



iV 2 (p)^',p) 



4 4 2.1 

-sinEa H — '-/aCosEcl sin2£'a 74Cos2£'a 

3 3 ' 3 3 ' . 



«(£>p) 



(33) 
(34) 



Notice that there are no finite-a errors here as required by current conservation. For the 
local current, it is found that 
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TABLE I. Simulation parameters. Njj is the number of gauge field configurations, a s t is the 
lattice spacing derived from the string tension, k c is the hopping parameter at the critical point 
and k s is the hopping parameter corresponding to the strange quark mass (from mx* /nix)- The 
error in the last digit is bracketed. 



Lattice 






a st [fm] k 


K c 










Improved Action 






6 3 xl2 


200 


6.25 


0.4 0.162, 0.168, 0.174 


0.1803(2) 


0.1658(8) 


8 3 xl4 


100 


6.8 


0.27 0.150, 0.154, 0.158 


0.1638(2) 


0.1561(5) 








Wilson Action 






6 3 xl2 


200 


4.5 


0.4 0.189, 0.201, 0.213 


0.2189(2) 


0.2050(3) 


8 3 xl4 


100 


5.5 


0.27 0.164, 0.172, 0.180 


0.1857(3) 


0.1779(7) 



(q(Z',p)\V 4 L \q(Z,p)) = 6& \l + 0(a 



(35) 



which again displays the desired classical improvement. 

Thus, although we are using an unimproved conserved current, the particular calculations 
being considered in this work are still improved to the level of the action. This conclusion will 
be valid for the numerical simulations as well, provided we create the external hadron states 
using operators that maintain the improvement. Local interpolating fields, such as those of 
Eqs. (|8|-|T0"D which will be used in the present work, guarantee that classical improvement is 
not destroyed. 



IV. NUMERICAL SIMULATIONS 

Calculations were performed at two values of (3 with the improved action, corresponding 
to lattice spacings of 0.4fm and 0.27fm as derived from the string tension. For comparison, 
the Wilson action was also used to calculate at the same lattice spacings. Some details of 
the simulations are provided in Table |. 

In all cases, we used pseudo-heatbath updating with periodic boundary conditions for 
the gauge fields. The lattices were thermalized with 4000 sweeps, and then 250 sweeps (200 
sweeps) were discarded between each pair of improved (Wilson) configurations that was kept. 
A stabilized biconjugate gradient algorithm was used for the fermion matrix inversion. To 
maximize the number of useful time slices, Dirichlet boundary conditions were used for 
fermions in the time direction, with periodic boundary conditions in all spacial directions. 
The source was placed two timesteps away from the time boundary in every case, and for 
the 3-point correlators of Rp, Ry and Rb we attached the vector current insertion to the 
lattice's centre time slice. 

Fig. |3| shows the effective mass plot for a pseudoscalar meson with the improved action. 
The non-monotonic behaviour near the source (t = 3 in Fig. is a consequence of next- 
nearest-neighbour interactions, as discussed by Luscher and Weisz [13]. The effect of these 



high-energy oscillations decreases as k increases, and Fig. shows that there is already an 
identifiable plateau at our smallest value of k. 
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FIG. 3. The effective mass plot for a pseudoscalar meson using the improved action at f3 = 6.8 
with k = 0.150. 



To extract the mass from an effective mass plot such as Fig. |3], or to obtain Zy or 1/fy 
from the appropriate plot, we require that the plateau contains at least three neighbouring 
points, and we never use the two time slices near the boundary (12 and 13 in Fig. |3|) 
which can be affected by the Dirichlet boundary condition. The systematic uncertainties 
associated with choosing a plateau region will not be presented since we do not expect these 
uncertainties to affect the principal conclusions of our study. All statistical errors will be 
computed using half the distance between the 16th and 84th percentiles from 500 bootstrap 
samples. 

The results for the renormalization factor, obtained from the improved action at (3 = 6.8, 
are shown for degenerate-mass quarks in Fig. [|. The value of Zy that comes from using 
meson decay amplitudes (Rd) m Eq. (|5]) is significantly smaller than the values that come 
from inserting a vector current between pion, p meson or nucleon states (Rp, Ry, Rb), 
and this difference is a reflection of the finite-a errors that are present in Eq. (||). To 
decide whether these finite-a errors should be considered as large or small, a comparison 
is made to results from the Wilson action. Fig. |] shows that the improved action results 
agree remarkably well with Wilson data at f3 = 6.0, while being clearly distinguishable from 
Wilson data (for R D ) at /3 = 5.85. 

By defining an effective hopping parameter which averages over the two quark flavours 
in a meson or baryon, 



2/t c ff 



2 V2ki + 2k 2 ) 



(36) 



we can plot all of our data, involving degenerate and nondegenerate quarks, together. Fig. || 
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FIG. 4. The renormalization constant as a function of the pseudoscalar mass squared with 
degenerate quarks. Solid symbols are from the improved action at = 6.8; open stars |14| and 



diamonds |15| are Wilson at (3 = 6.0; open squares [16| are Wilson at (3 = 5.85. 
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FIG. 5. The renormalization constant as a function of l/(2/i e ff). Open(solid) symbols are from 
the improved action at (3 = 6.25(6.8). 



11 



shows that our improved action results for Zy appear to be linear in l/(2« e fr)- 

The chiral limit, reached by extrapolating Zy, 1/fy, pseudoscalar meson 

masses squared, vector meson masses and baryon masses linearly in 1/k, using three data 
points (degenerate and nondegenerate quarks are used separately). In this way, we obtain 
the results of Fig. || for the renormalization factor without strange quark involvement. Also 
shown are the chiral limit values corresponding to the Wilson results of Fig. |], and our own 
Wilson calculations at smaller f3's. Note that the improved action approximately resembles 
the trajectory of Wilson data, but translated horizontally to lattice spacings about three 
times as large. This is very similar to what was also found in a comparison of the mass 
spectra J7|, and provides the first indication that the accuracy of calculations with tadpole- 
improved actions on coarse lattices persists beyond spectral quantities. 

Fig. |7| shows the p meson decay constant without including the effect of Zy, and with all 
Wilson quarks normalized by \^2k. Again, it is found that the improved action gives results 
which resemble the Wilson data, but at a lattice spacing about three times as large. 

Our best determination of the p meson decay constant comes from the product of Zy 
(obtained from Rp) and l/(Zyf p ), which is shown in Fig. ^. The Wilson calculations can 
be divided into two categories: those which compute Zy nonperturbatively as we have done 
in this work, and those which employ a perturbative estimate of Zy, denoted by Z pert . (The 
data in Fig. § with Z pert use quarks which are not normalized by the naive v2k, see the 
original papers for details.) The graph shows that both the improved and Wilson actions 
give nonperturbative results which are quite independent of lattice spacing over a large range 
of a, and which are near the experimental value. In contrast, the Wilson data using Z pert 
approaches a similar continuum limit but with a noticeable dependence on a. 

Having now discussed 1/fy in the chiral limit, let's return to our unextrapolated results, 
plotted in Fig. Notice that the slope of our data is consistent with the observed slope 
from the two experimental points l/f p and 1//^, but as a — > the simulations appear to 
favour values larger than the experimental ones. We find it interesting that the UKQCD 
data, also shown in Fig. || displays the same tendency. In their paper [19|, the UKQCD 



group concludes that this shift toward larger 1/fy values is the opposite of what would 
be expected from quenching, and they suggest that finite volume effects or lattice spacing 
effects may be responsible. The testing of such possibilities within the context of our 0(a 2 ) 
classically-improved, tadpole-improved action is much less computer-intensive than for the 
/? = 6.2 clover action, and might offer some valuable insight. 

The hopping parameter corresponding to the strange quark as given in Table |, was 
determined from the experimental relation 

rriK* , s 

— — = 1.8 37 
m K 

by linear interpolation in 1/k between the two neighbouring data points. Using this, Table [TI| 
gives our results for the three light vector meson decay constants and compares them to 
predictions from other groups. 
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FIG. 6. The renormalization constant in the chiral limit. Solid(open) symbols are from the 



improved(Wilson) action. The open star, diamond and square are from Refs. [jl4[|, [15| and [16] 
respectively. The lowest-lying point at each M p a is from the decay amplitude, Rd- 
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FIG. 7. The lattice decay constant in the chiral limit. Solid(open) symbols are from the im- 
proved (Wilson) action and the dashed line is the experimental value. Open squares, circles and 



diamonds are from Refs. |]lq| , [17| and [18| respectively. All Wilson quarks are normalized by V2k. 
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FIG. 8. The renormalized decay constant in the chiral limit. All symbols are defined as in 
Fig. [?| Squares and triangles use a nonperturbative determination of Zy, all other data rely on a 
perturbative estimate. 
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FIG. 9. Decay constant data before extrapolation/interpolation. Solid triangles (squares) are 
improved action data at (3 = 6.25 (6.8), open diamonds are (3 = 6.2 clover data [Q, and the 
asterisks represent experimental data. The horizontal scale is set by the physical p meson mass. 
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TABLE II. The decay constants from various simulations. Only the largest results from each 
reference are shown here. All clover and Wilson data use the perturbative renormalization constant 
Zpert- CL p is the lattice spacing derived from the p mass. 



Action 


a p [fm] 




!//*• 




iiiipiuvecL, jj — u.itj 
improved, = 6.8 
Wilson, = 6.0, Ref. @ 
Wilson, = 6.17, Ref. [17 
clover, = 6.2, Ref. fl§ 
Wilson, = 6.4, Ref. @ 
clover, = 6.4, Ref. [20| 
experiment 


0.242(5) 
0.0845(1) 
| 0.071(1) 

0.073(3) 

? 

? 


0.301(8) 
0.33(1) 
0.34(1) 

0.3161^ 
0.32(2) 
0.25(2) 
0.28 


0.284(5) 

0.298l| 
0.32(1) 
0.24(1) 


0.265(3) 

0.38(1) 
0.280lj? 

0.23 



V. CONCLUSIONS 

In this paper we have calculated vector current matrix elements and the renormalization 
factor in the light quark sector for an improved action on lattices with spacings of about 0.3 
to 0.4fm. This is the first step in going beyond the mass spectrum with improved actions 
on coarse lattices. 

The action used here is particularly simple. It has next-nearest-neighbour couplings to 
remove 0(a) and 0(a 2 ) errors at the classical level and also incorporates tadpole improve- 
ment. No special tuning of parameters and no field transformations are required with this 
action. 

Vector meson decay constants and the vector current renormalization factor calculated 
with the improved action were found to be comparable to those with the Wilson action on 
lattices with about 1/3 the spacing. This is the same as the improvement found previously 
in mass calculations and it represents a very substantial saving in computer resources. 
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